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I. INTRODUCTION 



In order to understand the dynamics of the very early universe, one must consider quan- 
tum gravity phenomena. With this aim, loop quantum cosmology (LQC) [l]-[3| confronts 
the quantization of symmetry reduced models in cosmology following loop quantum gravity 
(LQG) methods and ideas |4H6|. Although a full derivation of LQC from LQG has not yet 
been achieved, LQC is expected to correctly capture the behavior of the full theory, at least 
for those degrees of freedom which are responsible for the most important features of our 



universe. 



Most of the work to date in LQC has been devoted to the study of the so-called mini- 
superspace models, which only contain a finite number of degrees of freedom. Despite their 
simplicity, these models are surprisingly rich inasmuch as they cover many of the situations 
and phenomena of interest in LQC and have already shed light on Planck-scale physics in a 
cosmological setting. The study of the flat homogeneous and isotropic cosmological model 
in LQC [7 13] has shown that the classical big bang singularity is resolved by quantum 
geometry effects and that a quantum bounce occurs when the matter energy density reaches 
the critical density 0.41/?pi. All this happens in such a way that the large volume limit of 
the model (for semiclassical states) is well approximated by the classical dynamics dictated 
by general relativity, just as one should expect. In addition, further studies have shown that 
the covariant entropy bound holds in this model [jjj and that (assuming the presence of 
an inflaton field) the likelihood of obtaining viable initial conditions for slow roll inflation 



after the quantum bounce is higher than 99% [15] . It has been proven as well that the 



quantum bounce persists and that the large volume limit gives the correct classical limit 
when a nonzero cosmological constant is present [l6j and also for closed 17, 18|| and open 
[l~9| topologies. Finally, although less is known about anisotropic models as their quantum 
dynamics are considerably more complex, it has been shown that the classical singularity 
is resolved both in the spatially flat Bianchi type I model 2CH27| and the spatially curved 
Bianchi type II and type IX models 



One of the key features in LQC is that the curvature is expressed in terms of holonomies 
and one must determine how these holonomies are to be constructed. For homogeneous and 
isotropic space-times, it took several years to understand that one should use the so-called 
/2-scheme or improved dynamics scheme, presented in Ref. 0], rather than the /z -scheme or 
old dynamics scheme, which had been proposed in Ref. When anisotropic models were 
studied, there appeared two natural generalizations to the /2-scheme, which were presented 
in Ref. [23[ . Since one of these was considerably simpler to work with, and a priori seemed 
more natural than the other, a lot of the initial LQC work stud ying Bianchi I models 



26l . |32| . In particular, 



followed this procedure until several problems were pointed out 

the scaling properties of the (more complicated) alternate procedure have proven to be 
more suitable \2A 



27l |. Hence, unless one allows an explicit dependence of the /i-parameters 
on the coordinate cell adopted in the construction of the theory (or restricts the analysis 
exclusively to compact topologies where a distinguished choice of coordinate cell is available), 
the present consensus is that this alternate and more complicated procedure is the correct 
one to follow. We will refer to these two /2-schemes — i.e., the original simpler scheme and 
the more complicated one recently studied in Ref. (2?| — as schemes A and B respectively. 

The first aim of this paper is to complete the loop quantization of the vacuum Bianchi I 
model within scheme B, whose kinematical structure has been established in Ref. 27J and 
later analyzed in detail in Ref. [33[. We will see that the Hamiltonian constraint of the 
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model provides a difference equation in an internal discrete parameter v, which is strictly 
positive and proportional to the volume of the Bianchi I universe. Employing the form of 
the superselection sectors for the anisotropies, which were determined in Ref. j33j, we will 
show here that the quantum evolution equation is indeed well defined in vacuo, namely, that 
one can use v as a "time" variable and evolve the wave function in terms of it. In other 
words, we will show that a set of initial data, given on the section of minimum v, completely 
determines the physical solutions. Owing to this fact, we will be able to obtain the physical 
structure of the vacuum Bianchi I model for the first time. 

An extra motivation for the consideration of the vacuum Bianchi I model with the spatial 
topology of a three-torus is that its solutions coincide with the subset of homogeneous 
solutions (homogeneous sector) of the Gowdy T 3 model with linear polarization (34| . Based 
on the quantization of vacuum Bianchi I given here, one can then face the quantization of the 
Gowdy model in the framework of LQC, allowing for the introduction of inhomogeneities. 

The Gowdy T 3 model with linear polarization can be viewed as the simplest inhomo- 
geneous cosmological model. These cosmological spacetimes admit two axial Killing vector 
fields [lij and they describe universes devoid of matter which generically start with an initial 
curvature singularity [HBI, 36[ . Their quantization by standard methods has been discussed 



in detail in the literature (see, e.g., Refs. 37, 38j). It is now well known that after a complete 
deparametrization, the model admits an essentially unique Fock quantization with certain 
desired properties j38|. This explains the attention that has already been paid to the quan- 



tization of the Gowdy T 3 spacetimes in the framework of the loop theory [39|-|42[] . with the 
two-fold purpose of including inhomogeneities in LQC and achieving the quantum resolu- 
tion of the singularities of the model. In particular, Refs. [39|, |40j succeeded in doing this by 
proposing a hybrid quantization which combines the loop quantization of the Bianchi type 
I homogeneous sector -where scheme A was adopted for the improved dynamics- with the 
natural Fock quantization of Ref. [38[ for the inhomogeneities. Since it is generally accepted 
that scheme A must be replaced with the new scheme B in the quantization of the Bianchi 
I model, the second goal of this work is to discuss and complete the hybrid quantization of 
the Gowdy model making use of this alternate scheme for the homogeneous sector. With 
this aim, we will put on a rigorous basis some steps that are essential for the construction 
of a well-defined hybrid quantum model and which were left unfinished in Ref. [33| . 

Even though the kinematical Hilbert space of the hybrid procedure is naturally separated 
as the tensor product of two subspaces, namely the polymer space of the homogeneous sector 
times the Fock space for the inhomogeneities, the feasibility of this hybrid procedure is not 
at all trivial. As we will see, the Hamiltonian constraint of the model couples both sectors in 
a complicated way and it is not obvious whether the constraint can in fact be promoted to a 
well-defined operator, especially since the inhomogeneous sector contains an infinite number 
of degrees of freedom. Despite the complexity of the model, it was shown in Refs. 
that indeed the resulting hybrid quantization can be defined properly for scheme A and 
the corresponding physical Hilbert structure was obtained. Although a priori no relation 
of this structure with the kinematical one can be presumed, the standard quantum Fock 
description for the inhomogeneities was recovered in fact at the physical level. Now, the 
additional difficulties associated with scheme B, where holonomies along different directions 
no longer commute, make the new problem considerably more complicated. In this paper, 
we will prove that the scheme B hybrid quantization of the Gowdy T 3 model is viable. 

Let us mention that, apart from the analyses of the Gowdy model, other studies in the 
literature that have investigated the role of inhomogeneities in the framework of LQC are 
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Ref. (43[, which adopts a truncation of LQG and employs an approximation of the Born- 
Oppenheimer type, and the effective analysis of Ref. |44|. 

The paper is organized as follows. In Sec. [U we recall and summarize the kinematical 
structure of the vacuum Bianchi I model with T 3 topology in LQC adopting scheme B. In 
Sec. IHIt first we prove that the notion of evolution with respect to the volume is well posed as 
the associated initial value problem is well defined, and then we complete the quantization, 
characterizing the physical Hilbert space and a(n over) complete set of physical observables. 
In Sec. IIVI using the results obtained for the Bianchi I model, we show that the quantum 
dynamics of the Gowdy cosmology is also well defined when we employ scheme B in its 
hybrid quantization. Finally, we conclude in Sec. |V] with a discussion of our results and 
further comments. 



II. BIANCHI I T 3 MODEL IN VACUO: KINEMATICS 

This section summarizes the kinematical structure of the Bianchi I model quantized 
adopting scheme B for the improved dynamics prescription. We also include the description 



of the superselection sectors. We refer the reader to Refs. [27|, |33| for more details 



A. Vacuum Bianchi I Hamiltonian Constraint 

In order to describe the classical model, we choose angular coordinates 6, a, 5 G S 1 in 
which the spatial metric is diagonal. The elementary variables in loop quantum gravity are 
the Ashtekar-Barbero connection and the densitized triad. In the Bianchi I model, owing 
to homogeneity, each of them can be parametrized in a diagonal gauge by three coefficients. 
In terms of the fiducial co-triad {d9, da, d5} and the corresponding densitized fiducial triad, 
the coefficients of the Ashtekar-Barbero connection are given by Ci/(2n) and those of the 
densitized triad by pj/(47r 2 ), with % e {9, a, 6} and where the denominators of 2tt come from 
the periods of our coordinates. These coefficients form three pairs of canonical variables, 
with 

{c l ,p,} = 8nG 1 S lJ . (2.1) 

Here, 7 is the Immirzi parameter, G is Newton's constant, and Sij the Kronecker delta. 

For the study of Bianchi I models in LQC using scheme B, it proves useful to introduce 
the variables 



sgnjp^M and , 



(4tt7VA) 1 /3£ : 



pi 



pi 



PjPk 



(2.2) 



where £p\ = V Gh is the Planck length while the area A£ P1 is the gap in the spectrum of 
the area operator in LQG. In these equations, it is understood that repeated indices are not 
summed over and that the indices k are all different. Note that the above change of 
variables is in fact well defined only for nonvanishing triad variables Pi. Later on, we will 
see that this suffices for the study of the kinematical arena in LQC. 
It is also convenient to introduce the variable 

v = 2XeX a Xs, (2.3) 



which is proportional to the physical volume of the Bianchi I universe. 
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The operators that appear in the Hamiltonian constraint of the Bianchi I model for the 



considered scheme (see Ref. [27| for details) are A, and sinfej (or, equivalently, complex 
exponentials of It is easiest to work in the Aj representation, and in this case \ acts by 



multiplication while [45 



e^(A e ,A CT ,A 5 ) = # A 



Ao-, Xs 



(2.4) 



and likewise for the other complex exponentials of the fej's. Finally, it is most convenient to 
change the configuration variables from (A#, A CT , Xs) to e.g. (v,X a ,Xs), since v behaves in a 

simple manner under the action of e lbi , namely 



e ib °y{v, A (j, X s ) = * (v - 2 • sgn(A (T A 5 ), A CT , A a ) , 

&*{v, A (j, A tf ) = * f u - 2 • sgn^), - ~ 2 ' sgn(t;A - ) . A (j, A, 



(2.5) 
(2.6) 



where the symbol sgn denotes the sign function. The action of e lbs is obtained by inter- 
changing the roles of X a and Xs in the last equation. Thus, the effect of these holonomy 
operators on the dependence on v is just a constant shift (up to a sign). 

The Hamiltonian constraint operator obtained with a suitable choice of factor ordering 
and after a quantum densitization process 33| (such that the classical counterpart would 
appear in the action multiplied by the densitized lapse N = N/ \f\pePaPs\) is given in scheme 
Bby 

1 



c 



167rG7 2 



e CT + CT e + Q e Q s + 05©e + © CT ©5 + e 5 CT 



where 



sin(6i)sgn(Ai) + sgn(Aj)sin(6j 



(2.7) 



(2- 



Note that this is a different factor ordering choice for Gj than what was chosen in Ref. [2j 
This choice is preferable because the operator has the same action on wave functions sup- 
ported at large values of Aj while its action for small values of Aj is much simpler and the 
consequences are more transparent. Most importantly, the action of this operator does not 
allow any communication between different (Xg, X a , Xs) octants 25, 33, 39[. 

Under the action of Ch, the zero volume states decouple. Therefore, the singularities are 
resolved kinematically inasmuch as the quantum states that would correspond to them can 
be removed in practice from the kinematical Hilbert space Besides, under the 

action of this Hamiltonian constraint, the different octants remain invariant. Then, each 
octant contains different superselection sectors. This and the fact that the octants are all 
related by parity allow us to restrict our attention to the (strictly) positive octant. We do 
this for the remainder of this paper. Then, the Hamiltonian constraint acting on a wave 
function where v > 4 gives 



=v^ (v + 2)^h + lmt(v + 4, X a , X s ) -(v + 2)v^ X a , X s ) 



-{v- 2)v^^q (v, A CT , X s ) + (v - 2)V*>-4tf 4 (v - 4, A CT , A CT ) 



(2.9) 
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Here \E^ 4 are defined as follows: 

v + t x x s t ( v±n . v±2 \ f v±n 

± n, X a , X 5 ) = # [v ± n, —— ■ A CT , X s + * [v ± n, — — ■ A CT , A 5 

r f v±2 v±n \ f . v±2 x x 

+ * u ± n, A,j, — — • Aj + * v ± n, A CT , A, 

\ v v ±2 J \ v 

( v ±2 \ / i; ± n \ 

+ Vlv±n,X a , ■ X S J iv±n,X (T , ■ X 6 j . (2.10) 

On the other hand, if 2 < v < 4, the contribution of disappears in Eq. (12.9p . whereas 
if < f < 2 the two last contributions in that equation, namely those proportional to 
and \I/q , are absent. 



B. Superselection Sectors 



It was already pointed out in Ref. [27] that there are superselection sectors in v, denoted 
by a continuous parameter e G (0,4]. Given such a superselection sector, wave functions 
only have support on points where v — e + An, n being a natural number. 

Remarkably, there are also superselection sectors in the Aj's, although these sectors have 



a quite different structure compared to those in v. As it is shown in Ref. [33(, given e and 
an initial value A* > 0, the wave function will only have support on those points that can 
be expressed in the form 

*-(^)'n(S£)-* 

for some k, where mk, n^, and pk are nonnegative integers, and z is any integer unless e < 2, 
in which case z — 0. Note that different A*'s will yield the same superselection sector if they 
are related by Eq. (12. lip . It is not difficult to see that one of the available superselection 
sectors is the set of nonnegative rational numbers (it suffices to consider the case e = 2 with 
A* being a rational number). Furthermore, it follows from Eq. (12. lip that each superselection 
sector is countable, dense in M + , and that all superselection sectors are isomorphic 



Finally, one can show that the form of the superselection sectors for the three wave 
function variables (v, A CT , Xs) is characterized by three numbers (e, A*, A£), and is given by the 
tensor product of the superselection sectors for each individual variable, that is to say that 
there is no restriction on the sector of A^ given that of X a or vice versa [33j . In particular, 
note that if A* and A£ are compatible in the sense of Eq. (12.111) . then the superselection 
sectors of A CT and Xs are the same. It is clear that, given one superselection sector (e, A*, A£), 
a wave function will only have support on a countable number of points. Now, whereas the 
superselection sector in v only contains information about a discrete set of points separated 
by a constant shift, the superselection sector in the Aj's encodes the information of a set of 
points which are densely distributed in the positive quadrant. 



III. BIANCHI I T 3 MODEL IN VACUO: PHYSICAL STRUCTURE 



In this section we will analyze the solutions to the Hamiltonian constraint. We see from 
Eq. (12. 9p that the constraint provides a difference equation in the parameter v, and thus 
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it can be regarded as a (discrete) evolution equation in this parameter. In the previous 
analysis of Bianchi I carried out in Ref . 27J , a matter content was added, and the notion of 
evolution was developed in terms of a massless scalar field, instead of doing it in terms of this 
volume parameter v, which has a purely geometric nature. In the former case, thanks to the 
suitable properties of the massless field, which is quantized in a standard Schrodinger-like 
representation, it is straightforward to prove that the associated initial value problem is well 
posed. In fact, in this respect the situation is quite similar to that found in (relativistic) 
Quantum Mechanics. Nonetheless, regarding the geometry part, the physical structure of 
the solutions remained unanswered. Now, in the vacuum case considered here, the role of 
"time" is played by the volume variable v, which has been polymerically quantized. Because 
of its discrete nature, the fact that the associated notion of evolution is well defined is not 
trivial. We will show in this section that the dynamics are correctly posed: a set of initial 
data evaluated on the section of initial v completely determines the physical solution. As 
we will see, the proof is not direct, owing to the complexity of scheme B. In turn, this result 
will allow us to obtain for the first time the physical Hilbert space of the vacuum Bianchi I 
model in LQC and a(n over) complete set of observables, thus completing the quantization 
of the model. 



A. Solutions to the Hamiltonian Constraint 

Since we do not expect generic solutions to the Hamiltonian constraint to be normalizable 
in the kinematical Hilbert space, we will look for solutions in a larger space, namely the 
algebraic dual of a suitable domain of definition for the Hamiltonian constraint operator. 
It will be convenient to work with the variables Xj = InAj instead of the Aj's themselves, 
as the former variables run over the real line while the latter are positive and, besides, the 
Xi's suffer displacements under the action of the Hamiltonian constraint operator instead 
of dilatations or contractions. In fact, from Eq. (12. lip , the superselection sectors in Xi are 
formed by those points such that (x, — In A*) = Wi G Z e , where 

*-{. h (S^) + 5> ta (££)}. (3,) 

Here, for convenience, we have slightly changed the notation with respect to Eq. (12. lip , so 
that now rhk > fik are nonnegative integers, and pk can take any integer value. Recall that 
z is any integer unless e < 2 in which case z = 0. We note that Z e is dense in the real line, 
because the superselection sectors of Aj are dense in R + and the logarithm is a continuous 
function from the positive axis to the real line. In spite of the introduction of the Xj's, 
we will still keep v as one of our variables given its nice behavior under the action of the 
Hamiltonian constraint operator. 

Since the Wheeler-DeWitt equation associated to the Bianchi I model is actually a first 
order differential equation in the three variables x^, it should be possible to determine the 
entire solution to the Hamiltonian constraint supplying as initial data its restriction to 
one Cauchy slice (i.e., a surface with constant value of one of the Xj's or alternatively with 
constant v). However, the constraint in LQC is a second order difference operator for generic 
values of v. Therefore, it is not immediately clear how the solution can be determined from 
one slice of initial data obtained at a constant value of v. The solution lies in the different 
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form that the action of the Hamiltonian constraint operator has on states with v < 4, this 
form is that of a first order difference equation. 

Given a superselection sector in v, denoted by e G (0,4], one obtains this first order 
difference equation in v for v = e + 4 in terms of the initial data on the slice v — e. If 
one can solve this (highly coupled) difference equation, it will then be possible to solve the 
(again highly coupled) second order difference equation for v = e + 8 in terms of the data 
on the slices v — e and v = e + 4. One can then follow this strategy in order to obtain the 
full solution to the Hamiltonian constraint for all v. Finally, because the v + 4 terms always 
appear in the same combination, given by ^4 , we only need to show how to derive the data 
on v = e + 4, and then the data for all other v can be obtained in the same manner. 

The difference equation we are interested in is 



^2{e + 4:,x a ,xs) 



e + 4 



e - 2 

^0 (e, x a , x s ) + —7^0 ( e > x <n x s) 



(3.2) 



where the second term on the righthand side is absent if e < 2. Since the righthand side is 
known, the question is just whether one is able to obtain \I/(e+4, x a , x$) from ^4 (e+4, x a , x$) 
in order to derive the form of the wave function for all v. 
The explicit form of ^4 is given by: 

'v + 4" 



Vt(v + 4, x (j , x s ) =v(v + 4, x /j , In ^ 



+ ^ U + 4,x (T ,ln 



+ 2 
v + 2 



v 



+ x s ] - * I r - Lin 

+ x s 



+ * \ v + 4, In 
+ * ( v + 4, In 



v + 2 

v 

v + 4 
v + 2 



+ x u , In 
+ Xfj , In 



• M> ( u - - L In 

v + 4 



w + 4 
v + 2 
v + 2 



+ ^s 



+ x s 
+ x s 



(3.3) 



for any value of v (and, in particular, for v — e). It can be expressed as the result of the 
action of two separate operators, 



ty+(v + 4, Xfj , x s ) = U 6 (v + 4)AV(y, Xfj , x&), 
where A only shifts the value of v, namely, 

AV(v,x ,x s ) = (e 2 **} 2 V(v,x a ,x s ) = V(v + 4,x a ,x s ), 

while 



(3.4) 



(3.5) 



U 6 (v)^(v,x a ,xs) = 



—ibc 



+ e 



—ib a p—ib, 



■0 + e 



—iba p—ib, 



s + e -ib Se - 



+ e~ ibs e~ ibe + e~ ibe e~ ibs 



e ib e \ ^( Vi x a ,Xs) 



(3.6) 



has a trivial action on the ^-sector. 

The invertibility of the operator U 6 (v) for any value of v would guarantee that one 
can determine ^(v, x a , x$) from ^4 (v, x a , x$). Assuming for the moment that the inverse 
operator [Uq(v )] _1 exists, we can derive the state at the volume e + 4 by calculating 



^(e + 4,x a ,xs) 



e + 4 



-1 1 r e — 2 

£/ 6 (e + 4) ty+(e,x a ,x s ) + — — (e,x a ,x s ) 

L 6 -\- Z 



(3.7) 
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Once again, the second term in the square brackets on the right hand side of this equation 
does not appear if e < 2. 

Again assuming the existence of [^(f)] -1 , it is now straightforward to obtain the value 
of the wave function for the section v = e + 8: 



ty(6 + 8,X a ,Xs) 



e + 4 



U«(e + 8) 



e + 

e + 2 / < 



^+(e + 4, x a , x s ) + ^4*o (e + 4, x a , x & ) 
e + o 



e + 6 V e + 4 



^ 4 (e,x a ,xs) 



(3.8) 



It is clear how to repeat this procedure in order to get the value of the wave function for 
all larger v as well. Therefore, if the operator Uq{v) can be inverted we conclude that the 
initial value problem in terms of v is well posed, at least at v — e. In the next section we 
show that indeed this is the case. 



B. The operator Uq 

Let us then analyze the operator Uq(v) to see that its action can be inverted. First, 
we provide a suitable domain of definition for Uq(v), keeping fixed the value of v (in other 
words, we restrict the discussion just to a slice of constant v). For each direction i = o 
or 5, consider the linear span Cyl A * of the states whose support is just one point Xi of the 

i 

superselection sector determined by Eq. (13.11) . with (x; — InAj*) = Wi G Z e . We call H\* 
the Hilbert completion of this vector space with the discrete inner product. Then, we can 
choose the tensor product Cyl A * <g> Cyl A * as the domain for Uq{v). 
Now, if we define on the Hilbert space "H A * ® H\* the translations 

x (j , x s ) = *{v, w a + x (j , w s + x s ), (3.9) 

then the operator U§{v) is just a sum of six translations of this kind. These translation 
operators are unitary because, if w a and ws are two numbers in Z ei so that a shift of Xi 
by any of them leaves invariant the superselection sector [see Eq. (13.11) ]. then the sum of 
\U^ w,T ' Ws ' ) ^(v, x a , xs)\ 2 over all x a and xs in the superselection sector coincides with the sum 
of \ty(v , x a , xs)\ 2 . Moreover, owing to this property and the Schwarz inequality, we conclude 
that the norm of the operator t) 6 (v) is bounded by 6. 

Since Uq(v) is bounded, it can be extended as a well-defined operator to the entire Hilbert 
space. This extension [which we also denote by Uq(v)} provides in fact a normal operator 
— namely, the operator commutes with its adjoint — as the translations in x a and/or in xs 
commute. Hence, in particular, it is guaranteed that the residual spectrum is empty. Thus 
the operator U%{v) is invertible in our Hilbert space if and only if its point spectrum does 
not contain the zero. 

It is not difficult to convince oneself that the point spectrum of U$(v) must be empty 
owing to the properties of the operator. The idea is that since this operator is just a linear 
combination of translations, any of its eigenfunctions must possess a certain translational 
invariance which would prevent them from being normalizable. In order to see this, let us 
consider again the translations U( w '" w «\ with w a ,ws € Z t . Since they all commute with 
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each other as well as with U§{y), they can all be diagonalized simultaneously, that is, there 
exists a basis of common (generalized) eigenfunctions. Let us call in the following 

so that U^ Wc " WB ' = U™' T Up . Given one of the eigenfunctions common to all of these trans- 
lations, we denote the corresponding eigenvalue of Up by Pi(vJi) (with i = a, 5). This 
eigenvalue must be a complex number of unit norm, because the translation operators are 
unitary In addition, since Up Up = fjp +Wi , it follows that 

Pi(wi)pi(wi) = pi(wi + iBi). (3.11) 

Recalling that all points in the superselection sector can be reached from In A* by a trans- 
lation Up, it is a simple exercise to show that the eigenfunctions are proportional to 
Pa(w a )p5(ws). We can always change this wave function by a constant of unit norm, and 
thus we fix Pi(0) = 1. Besides, in order to determine completely the wave function, we only 
need to know the value of p%{wi) in an appropriate subset of Z ei namely any collection of 
noncommensurable points which can generate the entire set by multiplication by integers. 
It is possible to see that property (13. lip provides then all the information about pi at the 
rest of points in Z e . In particular, pi(nwi) = [pi(wi)] n . 

The wave functions Pi(wi) are clearly nonnormalizable with respect to the discrete inner 
product in "H A * as they have complex unit norm at each point of the superselection sector 
(the shift of Z e by In A*) and the sector contains an infinite number of points. In addition, 
different wave functions p«(wj) must be orthogonal, because there always exists a (unitary) 
translation operator on Tix* whose eigenvalue differs for the two wave functions. 

At this stage of the discussion, it is worth noticing that, by the very construction of 
the algebra of fundamental operators in LQC previous to the introduction of superselection 
sectors, the operators Up that act as translations in the %i representation can be identified 
in the holonomy/connection representation — where they act as multiplicative operators — 
as elements of the Bohr compactification of the real line, Msohr (HI- These elements can be 
understood as maps pi from the real line (corresponding to all possible real values of Xj, or 
equivalently of Wi) to the circle such that they satisfy condition (13.111) and Pi(0) = 1. Owing 
to superselection, however, the values of Wi are now restricted to belong to Z t . We can then 
identify the wave functions Pi(wi) as equivalence classes of elements in IReohr, the equivalence 
relation being the identification of all those maps pi which differ only by their action on the 
set complementary to Z^ in the real line, i.e., M\2 e . Examples of Pi(w{) are provided by the 
exponential maps exp (ikiWi) from Z E to S 1 . Since Z e contains noncommensurable numbers, 
these exponentials separate all real values of ki [that is, for any two values of ki one can find 
a value of for which the exponentials exp (ikiWi) are different]. So, the set of possible and 
distinct pi contains all the exponentials with ki G K. 

Returning to the operator U 6 (v), it is straightforward to find its eigenvalue for each of 
the analyzed wave functions. It is given by 

ui 6 (p„,Ps)=J2{Pt [ ln fe)]+p* M^i)]} 

i=<7,<5 

+ E (MW^)]P,[W^)]}- (3-12) 

i,j=a,6;i^j 
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Remember that here v > 4. The point spectrum of Ue(v) will not contain the zero pro- 
vided that there is no normalizable linear superposition of the above wave functions with 
cu6(Pa,Ps) = 0. In this superposition, the measure for pi is continuous: this is a consequence 
of the wave functions Pi{wj) not being normalizable in T-L\*. The restriction to the kernel of 

Ue(v) is achieved then by introducing a delta function of UQ(p a ,ps) (peaked at zero). If one 
computes the norm of this superposition, the orthogonality of the wave functions Pi(wi) leads 
to integrals over the square complex norm of each (p a , p^-contribution. But this contains 
a square delta, so that the norm diverges. Therefore, the point spectrum of the operator 
Uq(v) does not contain the zero, as we wanted to show. Actually, one can apply the same 
line of reasoning for any other possible eigenvalue of Uq, not just for zero, showing that in 
fact the point spectrum of this operator is empty. 



C. Physical Hilbert space 

Now that we have seen that the solutions to the Hamiltonian constraint are completely 
determined by the data on the initial slice v — e, we can identify these solutions with the 
corresponding initial data and characterize the physical Hilbert space by providing a Hilbert 
structure to the data, belonging in principle to the dual of the vector space Cyl A * ® Cyl A * C 

u K n K . ' ' ' ' 

In order to endow them with an inner product, we take a(n over) complete set of classical 
observables forming a closed algebra, and we impose that the quantum counterpart of their 
complex conjugation relations become adjointness relations between operators. Such a set is 
formed by the operators e lXi and Up, with G Z t and i = a, S. For i/j(\ a , Xs) € Cyl A * Cg>Cyl A * 
(and for the initial data by duality), these operators are defined as 

e^ip(x a , x s ) = e ix °ip(x a , x s ), (3.13) 
U^ip(x a , x s ) = ^(x a + to (j , x s ), (3.14) 

and similarly for e lXs and U^ 6 . Clearly, all these operators are unitary in , according 

with their reality conditions. Therefore, we conclude that this Hilbert space is precisely the 
physical Hilbert space of the vacuum Bianchi I model. 



IV. HYBRID QUANTIZATION OF THE GOWDY T 3 COSMOLOGIES 

The Gowdy T 3 model can be viewed as homogeneous Bianchi I backgrounds which allow 
certain inhomogeneous modes of the gravitational field to propagate along one direction. 
This natural separation in homogeneous and inhomogeneous sectors motivated a hybrid 
quantization of the model which combines the loop quantization of the Bianchi I phase 
space with a natural Fock quantization for the inhomogeneities, and which was carried out 



in Refs. |39|, |40| adopting scheme A for the improved dynamics in the quantization of the 
Bianchi I sector. This separation of degrees of freedom is nonperturbative and independent 
of the strength of the inhomogeneities at the classical level. Although ideally one should 
perform a LQC quantization for the inhomogeneous degrees of freedom as well, this hybrid 
approach is justified if the most relevant quantum geometry effects (but not necessarily all 
quantum effects) are those that affect the homogenous background so that one can establish 
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a kind of perturbative hierarchy in their treatment. In addition, it is natural to adopt a 
Fock quantization of the inhomogeneities in this context based on the expectation that a 
conventional Fock description of the inhomogeneities ought to be recovered from LQC in 
a regime where quantum geometry phenomena are negligible. In this case there exists a 
privileged Fock quantization under certain requirements on the symmetries of the vacuum 
and on the existence of a unitary dynamics with respect to an emergent time these 
properties provide a natural Fock quantization for the inhomogeneous modes. 

Here, we will show that the hybrid quantization of the Gowdy model employing scheme 
B for the loop quantization of the homogeneous sector is also viable. As in previous works 



with the other scheme [39l . I40j . this is not a trivial issue owing to the coupling between the 
homogeneous and inhomogeneous sectors in the Hamiltonian constraint. Moreover, now the 
structure of the homogeneous sector is much more complicated, owing to the intricacy of the 
holonomy operators in the new scheme B. The results obtained in Sec. II III will be essential 
in order to see that the hybrid quantization is well defined within scheme B as well. Our 



demonstration provides a necessary justification for the steps followed in Ref. [33| where the 
physical Hilbert space of this new hybrid Gowdy model was obtained. 



A. Kinematical structure and Hamiltonian constraint operator 



As in the Bianchi I model, since the spatial topology is that of a three-torus, we have 

one 



6, a, 5 G S 1 with a coordinate length of 2ir. Following a careful gauge- fixing [40J, |47 
finds that the information about the homogeneous degrees of freedom (which describe the 
subfamily of homogeneous space-times in the Gowdy model) can be encoded in the Bianchi 
I variables q and introduced in Sec. Ill Al On the other hand, the inhomogeneities cor- 
responding to the content of gravitational waves can be described by a single metric field 
(without a zero mode), which in turn can be described by creation and annihilation-like 
variables {(a m , a^J, m G Z — {0}}, defined in the same way as the natural variables that 
one would adopt if the field behaved as a free massless scalar field. Owing to the par- 
tial gauge-fixing, only two global constraints remain on the system: the zero mode of the 
Hamiltonian constraint, which generates time reparametrizations, and the zero mode of the 



#-diffeomorphism constraint, which generates translations around the ^-circle (see Ref. [40 
for details). 

In order to proceed with the hybrid quantization of the Gowdy model, we follow the 
LQC approach for the homogeneous degrees of freedom and, as in the Bianchi I model, we 
adopt the prescription pjCj — > ©; where the operators Oj are defined in Eq. (12. 8p while for 
the inhomogeneities we promote the creation and annihilation variables to operators in the 
standard quantum field theory fashion. The kinematical Hilbert space is then the tensor 
product of the polymer space of the Bianchi I model times the resulting Fock space for the 
inhomogeneities [33J. 

The generator of translations around the ^-circle only affects the inhomogeneities and it 



is straightforward to impose in the quantum theory [33l . l39l . |40|. On the other hand, as we 



pointed out earlier, the Hamiltonian constraint couples the homogeneous and inhomogeneous 



sectors in a nontrivial way. The resulting operator has the explicit form [33 
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Hf nt + 327r> e |H | , 



(4.1) 



where 



Ho = ^ Ha^a™ 

m^O 



and H L = J] 2|m[ ( 2a ™" m + 



Q'mfl'—m, ^Irfl—m 



(4.2) 



m^O 



are the terms that involve the inhomogeneities and the regulated "inverse triad" operator 
representing |pe| _1//4 is given by 33 



LW 1/4 J 



|f,A CT , As) 



'2 A ff Aa 



(4tt7VA^ 1 ) 1 /6 



v^k+i] - v>- i|) \v, a«t, a 5 ). 



(4.3) 



Note that the first term in Eq. f!4.2[) is the Hamiltonian of a free massless scalar field and 
the second is a quadratic interaction Hamiltonian. Owing to the coupling between these 
terms and those of the homogeneous sector, it is not guaranteed that the hybrid approach is 
physically feasible beyond the kinematical level, namely, once the constraints are imposed. 
In the remainder of this section we confirm that one attains in fact a well-defined physical 
theory. 

The explicit form of the action of the Hamiltonian constraint operator on kinematical 
states is easy to compute as most of the terms in the operator have already been considered 
in the Bianchi I model or are extensions of well-known operators. The most significant 

subtleties concern the action of HL. as this operator includes the sum of all a^aLm terms 



[m 



7^ 0), each of which creates an extra pair of "particles" in the modes m and — m. 

However, in spite of the fact that Hf nt creates an infinite number of particles, one can prove 
that it is a well-defined operator in a suitable dense domain of the Fock space |33| . On the 
other hand, it is straightforward to see that Ch leaves invariant Hilbert subspaces which 
are the tensor product of the superselection sectors of the Bianchi I model times the Fock 
space. Therefore, as in the Bianchi I model, we can restrict the study to separable Hilbert 
subspaces whose states have, in the homogeneous sector, quantum numbers (v, A CT , \$) with 
support in discrete sets contained in the positive octant. Let us remember that while v takes 
values in a semi-lattice of constant step equal to 4 with a minimum equal to e 6 (0,4], the 
values of (A CT , Xs) densely cover the positive quadrant of the real plane. 

Even though, from a physical perspective, one is only interested in small inhomogeneities 
which produce a perturbation around the homogeneous Bianchi I background, the hybrid 
quantum model is well defined and consistent without restrictions on the wave numbers or 
occupation numbers of the modes, and the evolution can be obtained in much the same 
manner as in the vacuum Bianchi I model. Using the result that the Bianchi I model in 
scheme B leads to a well posed initial value problem on the section of constant volume v — e 
from which one can evolve the physical state in steps of four units in v, one can show that 
the physical evolution in the hybrid Gowdy model is also (formally) solvable adopting a 
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perturbative approach, in which the effect of the interaction term H; nt is treated as small 

compared to the free Hamiltonian term Ho. Starting with initial data at v = e, one can then 
find the form of the physical wave functions at v = e + 4 in a perturbative expansion. With 
this data, one can continue the evolution to the next section v — e + 8. This procedure can 
be repeated until one obtains the expression of the physical wave function at the wanted 
value of v and up to the desired perturbative order. Actually, this perturbative expansion 
can be understood as an asymptotic expansion in the limit in which the Immirzi parameter 



tends to infinity. The details of this perturbative expansion are presented in Ref. [33 . 

The important point here is that the evolution is well defined in this perturbative ap- 
proach. This is mainly due to the fact that the initial value problem for the vacuum Bianchi 
I sector is well posed and therefore this result depends upon the proof presented in Sec. III. 
In this sense, the initial value problems in the vacuum Bianchi I model and the Gowdy T 3 
model are closely related. 



V. DISCUSSION 

In this paper we have first considered vacuum Bianchi I universes with a three-torus 
topology in the framework of LQC, adopting a new scheme for the improved dynamics which 



was put forward in Ref. [271 ] . We have examined some of the aspects of this quantization 
which had remained unanswered in Ref. (27[, like the decoupling of triad components with 
different orientations under the action of the Hamiltonian constraint, the structure of the 
superselection sectors in the anisotropies, and the evolution of physical states in terms of 
the volume as a discrete, internal evolution variable. Then, we have shown that the initial 
value problem is well posed and completed the quantization of the vacuum Bianchi I model, 
following scheme B for the implementation of the improved dynamics. 

In Sec. IIVI we have used the results regarding the vacuum Bianchi I model in order to 
show that the scheme B hybrid quantization of the linearly polarized Gowdy T 3 cosmological 
model is viable. This hybrid quantization provides a first step towards a better understand- 
ing of the effect of inhomogeneities in LQC; this is necessary if one wants to eventually 
obtain predictions about the influence and possible traces of quantum gravity in phenomena 
like primordial gravitational waves, the cosmic microwave background, and the physics of 
the early universe in general. 

The loop quantization of the Bianchi I model leads to superselection in separable sectors 
not only for the volume, but also for the anisotropies. Moreover, every superselection sector 
is restricted to an octant. This is because the Hamiltonian constraint operator, due to appro- 
priate factor ordering choices, does not mix eigenstates of the densitized triad components 
with different orientations. Moreover, while the superselection sectors in the volume of the 
Bianchi I universes consist of equidistant points forming a semi-lattice, the superselection 
sectors in the anisotropies are dense sets in the real semi-axis. On the other hand, the 
restriction to a definite orientation of the triad components without imposing any kind of 
boundary conditions, together with the fact that the initial value problem for the evolution 
is well posed at the minimum value of v (i.e., v = e), can be regarded as a realization of a 
no-boundary prescription for the dynamics. In addition, it is worth emphasizing the result 
that the discrete evolution in v is well defined starting from the initial section v — e. If this 
were not the case, the evolution would break down for Bianchi I cosmologies in vacuo for 
scheme B and, without any reasonable justification, the inclusion of matter would turn out 
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to be critical in order for the dynamics of the model to be viable; we have shown that this 
is not the case. 

Concerning the hybrid quantization of the Gowdy model, an important point is that 
the LQC/Fock split that we have considered assumes that the quantum behavior of the 
inhomogeneities can be well approximated by conventional quantum field theory methods 
so that any quantum geometry effects due to the presence of these inhomogeneities can 
be neglected as perturbatively small. In a true loop quantization of all the gravitational 
degrees of freedom (i.e., presumably in a reduction of LQG by a suitable incorporation 
of the symmetries of the Gowdy model), one should treat both the homogeneous and the 
inhomogeneous sectors in the same manner, that is to say, all of the degrees of freedom 
should be quantized a la loop. From this perspective, the main assumption in our analysis 
is that the qualitative results of the hybrid quantization capture the physics of a full loop 
quantization so long as the inhomogeneities are not directly affected in a significant way by 
quantum geometry phenomena. It is worth noticing that similar assumptions are implicit 
in the treatment of other models with matter in LQC inasmuch as matter fields are usually 
quantized by standard methods rather than by adopting a unified polymer quantization for 
all of the degrees of freedom, gravitational or not. 

Notice nonetheless that, even at this level, one can see that the hybrid quantization 
approach is sufficient to ensure that the classical cosmological singularities are resolved as 
the singular states corresponding to vanishing Bianchi I scale factors decouple under the 



action of the Hamiltonian constraint operator [33l. |48| . In addition, the quantum dynamics 
are well posed as one can use v as an evolution variable and, given the wave function at 
v = e, one can derive the wave function for all other values of v in the same superselection 
sector in a perturbative expansion in the interaction term for the inhomogeneities. 

There remain many open questions to be addressed, of course, the most important being 
a numerical study of the evolution of the wave function. This is a very difficult task as even 
the vacuum Bianchi I model in LQC has not yet been studied numerically in scheme B. We 
propose to begin with a simpler task and study effective equations associated to the model; 
this should yield some insight into the most relevant quantum geometry corrections to the 
classical model [49[. The most important point, however, is to understand inhomogeneities 
in LQC more deeply and to do this one will have to consider more general inhomogeneous 
space-times within the framework of LQC in order to obtain physical predictions about our 
early universe and understand their consequences. 
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